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in 
o 

The coalescence rate of two neutron stars (NS) is revisited. For estimation of the number 
^ ■ of bound NS-NS and the probability of their coalescence in a timescale r, the galactic star 

' formation history, directly derived from observations, and the evolution of massive stars 

^ I are considered. The newly established galactic merging rate is (1.7 ± 1.0) X 10^^ V"^^ , 

, while the local merging rate, including the contribution of elliptical galaxies, is about a 

factor of two higher, 3.4 X 10~^ yi'~^ ■ Using the present data basis on galaxy distribution 
in the local universe and the expected sensitivity of the first generation of laser beam 
interferometers, we estimate that one event should occur every 125 years for LIGO 
, and one event each 148 years for VIRGO. The situation is considerably improved for 

■ advanced-LIGO since we predict that 6 events per year should be detected whereas for 

a recently proposed VIRGO new configuration, the event rate might increase up to 3 
events every two years. 

Keywords: gravitational waves; neutron star binaries; laser interferometers 



1. Introduction 

The merger of two neutron stars (NS-NS), two black holes (BH-BH ) or a black 
hole and a neutron star (BH-NS) are among the most important sources of grav- 
itational waves, due to the huge energy released in the process. In particular, the 
coalescence of NS-NS binaries may radiate about 10^"^ erg in the last seconds of 
their inspiral trajectory, at frequencies up to 1.4-1.6 kHz, range covered by most 
of the ground-based laser interferometers like VIRGO\ LIGO^, GEO^ or TAMA". 
For understanding of the process dynamics, the prediction of its signal waveform, 
the estimate of its released power and the event rate, theoretical studies on the 
coalescence of NS-NS systems have been performed in the past years^^^'^. The dy- 
namics of the inspiral and, in particular, the study of tidal effects in the very last 
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orbits, may put strong constraints on the equation of state of dense matter. Con- 
cerning the coalescence rate of NS-NS binaries, theoretical evaluations are always 
performed in two main steps: firstly, the merging rate in our Galaxy is estimated 
and then, assuming that this rate is typical, an estimate is performed for the vol- 
ume of the universe sampled by a given detector, using some adequate scaling. The 
galactic merging rate has been estimated either by using binary population synthe- 
sis modcls^'^'"'^^ or from the statistics of the observed NS-NS binaries"'^'^'"'^^'^'''^^'"'^^. 
These estimates may differ by one order of magnitude, ranging in general from 10~^ 
yr~^ up to few times 10~^ yr~^, but values as high as ~ 3 x 10~'^ yr~^ have been 
reported in the literature^^'^^ . Recently Ref. 18 included the new highly relativistic 
NS-NS binary system PSR J0737-3039 in their statistical analysis, obtaining a value 
(~ 8.3 X 10~^ yr~^) comparable to the high rates derived previously by Ref. 16 and 
Ref. 17. Once the present galactic merging rate has been evaluated, the expected 
detection rate can be estimated by scaling the total luminosity within the volume 
probed by the detector with respect to the luminosity of the Galaxy at the same 
wave-band. Within the local universe (z < 0.01), the distribution of galaxies is not 
homogeneous and if the actual distribution is not taken into account, the expected 
detection rate may be seriously underestimated'^^. For instance, the inclusion of the 
huge concentration of galaxies in the direction of Norma-Centaurus may consider- 
ably increase the detection rate of these events^°. Estimates of the galactic merging 
rate performed according to a given evolutionary model, require the knowledge 
of the star formation history. Thc;sc> models assume visually that the star formation 
rate is proportional to the available mass of gas and a consequent exponential decay 
with time. However, observations and numerical simulations suggest that "disky" 
galaxies, including the Milky Way, are gradually formed by accretion and merger 
episodes^ ''^^, which affect the star formation history. It is clear that in spite of all 
efforts made in the past years to evaluate the local coalescence rate and the ex- 
pected detection rate by a given interferometric detector, the uncertainties are still 
quite large with consequences bearing on the future development of instruments and 
the detection strategy. In the present work, we combine the advantages of different 
approaches to present a new estimate of galactic merging rate based on population 
synthesis of the pulsar population^'''^^ and on new simulations on the; evolution 
of massive binaries. From the latter, we have derived the distribution probability 
P(t), which gives the fraction of newly formed NS-NS binaries with a coalescence 
timescale r, due to gravitational radiation losses. For the first time the star forma- 
tion history of our Galaxy, derived directly from observations, was included in the 
computation of the galactic NS-NS merging rate. Since elliptical galaxies contribute 
also to the morphological composition of galaxies in the local universe, we have esti- 
mated the expected coalescence rate in these objects, by adopting a star formation 
model able to reproduce their observed photometric properties. Then, the local av- 
erage coalescence rate has been estimated and weighted according to the total light 
fraction contribution of each morphological type. In a second step, we evaluate the 
volume of the universe probed by a given detector and compute the total luminosity 
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inside this volume, using available data basis. The organization of this paper is the 
following: in Section II we present the derivation of the local coalescence rate, in 
Section III wc estimate the expected detection rates taking into account the planned 
sensibility of the different interferometers and, finally, in Section IV we discuss our 
results and summarize our main conclusions. 

2. The local coalescence rate 
2.1. The galactic rate 

In order to calculate the coalescence rate of NS-NS pairs in a given galaxy, we 
adopt here the approach by Ref. 20. Let us suppose that a massive binary (masses 
of components higher than 9 Mq) is formed at the instant t' . Let r* be the mean 
evolutionary timescale required for the system to evolve into two neutron stars, 
typically of the order of lO'' — 10^ yr. Define P{t) as the probability per unit of 
time for a newly formed NS-NS binary to coalesce in a timescale t and define 
R^,{t) as the star formation rate, given in M0.yr~^. Under these conditions, the 
coalescence rate at instant t is 



where fb is the fraction of massive binaries formed among all stars, /3„s is the fraction 
of formed binaries which remain bounded after the second supernova event and A is 
the fraction per unit mass of formed stars in the mass interval 9-40 Mq . We assume 
that progenitors with initial masses above 40 Mq will produce black holes. We have 
set To as the minimum timescale for a NS-NS binary to coalesce. This minimum 
timescale was estimated from simulations to be described below. If the initial mass 
function (IMF) is of the form ^(m) = Am~^ , with 7 w 2.35 (Salpeter's law), 
normalized within the mass interval 0.1 - 80 Mq such as J m^{m)dm = 1, then 
A = J* 0£_{m)dm = 5.72 x 10""^ -^0^- The parameters (3ns, tq and the probability 
distribution P{t) were calculated from numerical simulations described in detail 
by Ref. 25. Here we outline only the main aspects of these simulations. A massive 
binary is initially generated according to the following prescriptions. The mass of 
the primary is obtained from a probability distribution corresponding to a Salpeter's 
IMF, while the secondary mass is derived from the observed mass ratio distribution 
for massive binaries obtained from observations The initial pair separation a 
is fixed by a probability distribution P (a) cc da /a '^^ , normalized between Rmm 
and Rmaa;- The minimum separation was taken to be equal to twice the Roche lobe 
of the primary and the maximum separation was taken to be equal to lOORmin- 
The initial eccentricity of the orbit was assumed to obey a distribution probability 
P{e)de = 2ede, corresponding to orbits filling completely the phase space. Until 
the explosion of the first supernova (the primary star), the orbital parameters of 
the system vary due to the mass-loss produced by the stellar wind from both stars 
or mass transfer. The latter mechanism is more rare, since these massive stars do 
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not frequently reach or overlap the Roche-lobe. The explosion of the first supernova 
leaves a NS remnant of 1.4 Mq and produces an "instantaneous" mass loss, which 
may disrupt (total orbital energy positive after the event) or not (total orbital 
energy negative after the explosion) the originally bound pair. Our simulations 
indicate that 77.4% of the systems remain bounded after the first supernova. The 
bound binary, now constituted by a NS star and an evolved star, which has already 
lost part of its original mass, enters into a second phase of slow mass-loss, with the 
orbital parameters varying as before, until the second supernova explosion occurs. 
We recall that our evolutionary scenario is similar to that developed by Ref. 10, 
in which none of the stars ever had the chance of being recycled by accretion. 
We have performed 500,000 numerical experiments, from which resulted 11,627 
NS-NS bound pairs, implying P^s = 0.024. A velocity kick is probably imparted 
to the nascent neutron star although its mechanism is still a matter of debate. 
Natal kicks may unbind binaries which otherwise might have remained bound or, 
less probably, conserve bound systems which without the kick would have been 
disrupted. In these experiments we have assumed that neutron stars have a natal 
kick velocity corresponding to a 1-D velocity dispersion of about 80 km/s. We have 
also investigated how these results are changed if a Maxwellian velocity distribution 
with a dispersion 1-D of 230 km/s is adopted. In this case, the resulting fraction 
of bound binaries is reduced by one order of magnitude, e.g., /3„s = 0.0029. These 
numbers corresponding to the resulting fraction of bound systems are consistent 
with previous analyses on effects of the natal kick Recent investigations on the 
spin period-eccentricity relation for NS-NS systems suggest that such a correlation 
can only be obtained if the second neutron star receives a kick substantially smaller 
(velocity dispersion less than 50 km/s) than kick velocities commonly assumed 
for single radio pulsars. If this conclusion is correct, values of (3ns derived from 
simulations with low kick velocity dispersion are more likely and will be adopted 
here. Clearly, the natal kick amplitude remains the major source of uncertainty in 
the estimate of the fraction of bound NS-NS binaries. However, once the relative 
number of NS-NS binaries is fixed, the fraction /(, of massive binaries formed among 
all stars is related to /3„s through the equation 



where Np and Nt, are respectively the number of single pulsars and the number of 
NS-NS binaries in the Galaxy. We have taken into account that the population of 
single pulsars results not only from the evolution of single stars but also partially 
from disrupted massive binaries. Since /h = 1 is a strong upper limit, this imposes a 
lower limit for (3ns for a given ratio Np/Nu. The number of isolated pulsars Np has 
been derived from population synthesis of single radio pulsars. We have employed 
the code described in detail in Ref. 23, 24, up-graded to take into account the new 
pulsars discovered at high frequencies by the Parkes Multibeam Survey as well as 
the most recent model for the pulsar velocity distribution The method consists 
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in generating population of pulsars with a given set of birth properties, following 
their evolution according to the dipole braking model and modeling selection effects 
that constrain radio detection. The best agreement between simulated and observed 
distributions of physical properties as the period, period derivative, distance, among 
others fixes the initial period and magnetic braking timescale distributions (see 
parameters in Table 1), as well as the total number of pulsars in the Galaxy {Np » 
250000) and their birthrate (one pulsar every 90 yr). Simulations using the recent 
results by Ref. 31 for the velocity distribution of radio pulsars were also performed 
and no significant differences in the values derived for parameters shown in Table 1 
or for the total number of pulsars were noticed. 

Table 1. parameters of the initial period 
Pq and magnetic braking timescale In to 
distribution, assumed to be Gaussian 

mean dispersion 

Po(ms) = 240±20 o-po= 80 ±20 
Inrois) = 11 ± 0.5 ainro = 3-6 ± 0.2 



The number of NS-NS binaries can be derived by linking the properties of the 
youngest pulsar in a given double neutron star system, supposed not to have been 
affected by external torques, to those of the population of single pulsars. In other 
words, if the evolution of the newly formed (second) pulsar was not yet affected by 
torques other than the canonical magnetic dipole, its rotation period increases as 
those of isolated pulsars. For this set of simulations, we have adopted a velocity dis- 
tribution for the baryccnter of the NS-NS binaries different from the velocity distri- 
bution adopted for single pulsars. We consider that millisecond pulsars (MSP) were 
spun-up by mass transfer, before the second supernova explosion, which disrupted 
the system. According to Ref. 32 MSP have a Maxwellian velocity distribution, cor- 
responding to a mean space velocity of 130 km/s. This velocity distribution must be 
that of those binary systems that survived the first supernova explosion There- 
fore, wc have supposed that binary systems which have also survived to the second 
explosion have a velocity distribution similar to those which have survived to the 
first explosion. The mean transverse velocity of observed MSP is around 90 km/s, 
which corresponds exactly to the Maxwellian mean space velocity of 130 km/s de- 
rived by Ref. 29 from parametric likelihood analysis techniques. So far, two binary 
systems have been detected through the second pulsar: PSR B2303-I-46 and PSR 
J0737-3039 B. In order to observe two pulsars having millisecond like kick velocities, 
our simulations predict a total number of binaries Nh = 730, corresponding to a 
ratio Np/Nb = 342. The advantage of our method is that we don't have to infer the 
statistical properties of binary pulsars from the very small observed sample and the 
only data we care about is the number of observed second born pulsars. The derived 
ratio Np/Nb implies from Eq. (2), ft, = 0.136 if the value /3„s = 0.024 is adopted. 
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Notice that if the upper hmit = 1 is imposed to Eq. (2), then we conclude that 
Pns > 0.0058 or that the 3-D velocity dispersion of the kick distribution of binaries 
should satisfy a < 270 km/s, e.g., lower than that of the single population. 

For each simulated NS-NS binary, the separation, the orbital period and the or- 
bital eccentricity are available, allowing the calculation of the coalescence timescale 
T due to emission of gravitational radiation. The cumulative number of pairs with a 
coalescence timescale less than r varies as N{< r) oc lg{T), corresponding to a prob- 
ability distribution P(r) = B/t (see Fig. 1) The simulations indicate a minimum 
coalescence timescale tq = 2x10^ yr and a considerable number of systems having 
coalescence timescales higher than the Hubble time. The normalized probability 
P(r) in the range 2 x 10'"^ up to 20 Gyr impHes B = 0.087. 

The other important term required for the evaluation of the coalescence rate from 
Eq.(l) is the star formation rate R*(f). In most of population synthesis models, in 
the absence of a detailed star formation theory, (t) is assumed to be proportional 
to the available mass of gas and under this (usual) assumption, from the mass- 
balance equation, it results i?*(i) oc e~"*. In fact, we will use such an approach to 
model the evolution of an elliptical galaxy, as we shall see in the next section. In 
some even more simplified galactic models the star formation rate is assumed to 
be constant. However, for the Galaxy, the star formation history can presently be 
derived from observations and most of recent studies have inferred that the star 
formation activity is non-monotonic with time (Ref. 35 and references therein) . Re- 
cent studies on the chromospheric activity index versus age relationship permitted 
the determination of ages for a sample of 552 stars from which was possible to re- 
construct the star formation history of our Galaxy^^. These data indicate enhanced 
periods of star formation at 1 Gyr, 2-5 Gyr and 7-9 Gyr ago, probably associated 
with accretion and merger episodes from which the disk grows and acquires angu- 
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Fig. 1. The cumulative fraction of simulated NS-NS binary systems with coalescence time less or 
equal t. Time is given is years and the solid curve shows the best fit to the simulated data. 
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lar momentum^^ . Here we use the results by Ref. 35 to compute numerically the 
integral defined in Eq. (1). Using the numbers obtained, it results for the present 
galactic NS-NS coalescence rate vs = (1.7 ± 1.0) x 10~^ yr~^. The estimated error 
is mostly due to uncertainties in the ratio Np/Nh derived from simulations. Simi- 
lar population synthesis calculations to estimate the birth-rate of compact binaries 
have also been performed by Ref. 36. However some differences exist between their 
and the present calculations. In the one hand, Ref. 36 normalized their birthrates to 
the type Il-f-Ib/c supernova frequency, which was assumed to be constant over the 
lifetime of the galactic disk. Moreover, they have assumed that NS-NS binaries have 
the same natal kick distribution as that observed for single pulsars. On the other 
hand, our calculations took into account the star formation history of our Galaxy 
and the normalization was derived from two simulation sets, which fixed respectively 
the fraction of bound NS-NS pairs and the ratio between binary and single pulsars. 
However, it is worth mentioning that the cumulative number of NS-NS pairs with a 
coalescence timescale less than r, resulting from both set of simulations is similar. 

2.2. The coalescence rate in ellipticals 

Excluding the Milk Way, no direct information is available for the neutron star pop- 
ulation and the star formation history of other galaxies. Thus, the coalescence rate 
in extragalactic objects can only be estimated theoretically. In this case, we have 
assumed for elliptical galaxies the same values derived for the Milky Way for the 
parameters fj,, fins- It is worth mentioning that previous estimates of the present co- 
alescent rate in ellipticals^^ lead to values which were about a half of that estimated 
for the Galaxy. However, in that work the star formation efficiency was constrained 
essentially by the present amount of gas in ellipticals, a very uncertainty quantity. 
A grid of models for elliptical galaxies was recently built in Ref. 37, including the 
effects of mass-loss by a galactic wind, responsible for the enrichment of the intra- 
cluster medium in metals. In these models, the parameters characterizing the star 
formation efficency and the IMF were chosen by an iterative procedure in order to 
reproduce the color-magnitude diagram, (U-V) against My of ellipticals in Coma 
and Virgo clusters. We take their model 3 as a representative of a typical E-galaxy, 
defined by an initial mass equal to 2x10^^ M© and a present luminosity Lb = 
2.9x10^° Lb 0. The IMF in these models are slightly flatter than the Salpeter's law 
(7 f« 2.19) and the corresponding A parameter is A = 8.71x10"^ ^q^- It should 
be emphasized that although the color-magnitude diagram imposes essentially con- 
straints on low mass stars, the models also explain quite well the relative abundance 
of the a-elements, determined by massive stars, since they reproduce the observed 
metallicity indices like Mg2 and < Fe > Under these conditions, the estimated 
present coalescence rate is ve = 8.6x10"^ yr~^. Notice that this value is about a 
factor of 5 higher than the galactic rate estimated in the previous section. Thus, in 
spite of fact that today ellipticals have practically a non-existent (or a very low) 
star formation activity, since the bulk of their stars were formed in the first 1-2 
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Gyr, the pairs merging today were copiously formed in the past with long coales- 
cence timescales. A similar phenomenon occurs with type la supernovae. No type II 
supcrnovac arc detected in ellipticals because these objects are the consequence of 
the evolution of massive stars, which have a rather short lifetime. However, type la 
supernovae are the only class detected in early type galaxies. This class of supernova 
is the consequence of the evolution of binary systems constituted by intermediate 
mass progenitors, which evolve into white dwarfs. The explosion of the compact 
object depends on the distribution of the merging timescale and only systems with 
long evolutionary timescales will explose by now. 

In Fig. 2 we show the evolution of the merging rate for the Milky Way and for our 
representative model of an elliptical galaxy. The smoothness observed in the curve 
describing the evolution of the coalescence rate in the E galaxy is a consequence 
of the fact that the bulk of the stars are quite old, being formed in a short time 
interval. Since the probability to have a merger is inversely proportional to the time, 
the curve reflects essentially this fact. This is not the case for our Galaxy, since its 
disk is built gradually and the star formation rate may suddenly vary as new matter 
is accreted. 



2.3. The local coalescence rate 

The models developed by Rcf. 37 can equally reproduce the colors of SO galaxies. 
If we assume that the coalescence rate estimated above is typical for E and SO 
galaxies having absolute magnitudes of about = -20.7 and that the fraction of 
these objects is about 35%, then the weighted local coalescence rate can be written 




Fig. 2. Evolution of coalescence rates for the Milk Way and for a typical E galaxy. Notice that 
in the elliptical the bulk of stars is formed very early and the merging rate reflects mainly the 
probability of a NS-NS binary to coalesce in a given timescale. For the Galaxy, the star formation 
is continuous and intermittent, producing a modulation effect in the coalescence rate history 
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as 

l^c = Mfs + fE'^^) (3) 

where fs = 0.65 and /b = 0.35 are respectively the adopted fractions of spirals 
and (E+SO) galaxies, the ViS and Lj's are the respective coalescence rates and lu- 
minosities. ^From this relation and values obtained before, the local mean weighted 
coalescence rate is Uc = 3.4x 10~^ yr~^. In the one hand, had we usually assumed the 
local merger rate equal to the galactic one, neglecting the contribution of (E+SO) 
galaxies, we would have obtained a value smaller by a factor of two. On the other 
hand, the very recent galactic rate derived by Ref. 18 seems to overestimate by 
a factor of 2.5 our results for the total local rate, which include NS-NS mergers 
occurring in galaxies other than the Milky Way. 

3. Expected detection rates 

Extrapolation from the local coalescence rate to the expected rate within the volume 
is made multiplying the local value by the factor Kb{< D), which is defined 
as the ratio between the total blue luminosity within the considered volume and 
the Milky Way luminosity {Lb,mw ~ 2.3 x 10^° L©). If galaxies were distributed 
homogeneously in the local universe, the scale factor correcting the local coalescence 
rate would simply be Kb{< D) k, O.Q2D\f^^. This approximation is probably rea- 
sonable for distances beyond 150 Mpc but underestimates the total luminosity in 
the B-band at lower distances, due to significant anisotropies in the distribution of 
galaxies. Moreover, the sensitivity of interferometric detectors is not isotropic, being 
larger to signals coming from above and below the plane. In this paper, the direc- 
tional sensitivity of the detector and the anisotropy in the distribution of galaxies 
will not be taken into account as in Ref. 19. However, real counts of galaxies inside 
a given volume will be considered, allowing a more realistic evaluation of the total 
luminosity. Analyses of the spatial distribution of galaxies for z < 0.033 have been 
performed by Ref. 38 (and references therein) using the LEDA (Lyon-Meudon Ex- 
tragalactic Database) catalog, which contains over one million of galaxies covering 
all-sky, including also a sub-sample of 134.000 objects having a measiircd rcdshift. 
The completeness of the catalog was evaluated by counts of galaxies up to a given 
magnitude, excluding galaxies in the zone of avoidance (b<| 15° |), which includes 
superclusters like the Great Attractor. According to Ref. 38, the redshift catalog 
is complete at level of 84% up to B = 14.5 and at a level of 90% for galaxies 
brighter than B = 11.0. It is worth mentioning that up to B = 14.5, the total lu- 
minosity does not increase as the cube of the distance, but more slowly with an 
exponent ^ 2.5. Ref. 38 considered different volume limited samples, including only 
galaxies brighter than a given absolute magnitude. We correct for the absence of 
faint galaxies, supposing a Schechter luminosity function in the magnitude interval 
— 15 > Mb > —22. Recent deep optical surveys in the zone of avoidance ^^'^O'^i 
confirm the existence of a large structure dubbed the "Great Attractor " in the 
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region | 6 |< 10° and 316° < I < 338°. This structure seems to be the intersection 
between the Centaurus Wall and the Norma supercluster, including some other 
clusters as the Centaurus-Crux complex. The center of the Great Attractor seems 
to be located at the Norma cluster (A3627) at Vz = 4844 km/s "^^ and up to now 
4423 galaxies have been identified as members of this supercluster. The total mass 
of such a complex is uncertain. A lower limit of 4x10^^ Mq is derived from direct 
counts and the observed velocity dispersion of the clusters, under the assumption 
of dynamical equilibrium. From large scale motions in this region, a mass as high 
as SxIO^^Mq was obtained^^. If we assume a more conservative value of about 
10^^ Mq for the total mass of the complex and a typical mass-to-luminosity ratio of 
150 Mq/Lb,q, then the Great Attractor alone would give a relative contribution to 
the total luminosity of Kb 2900, at a distance of 70 Mpc. This contribution was 
added to the total luminosity derived from data for distances D > 70 Mpc. In Fig. 
3, the calculated NS-NS coalescence rate ratio inside a volume of radius D, given 
by local rate times the factor Kb{< D), is shown as a function of the distance D. 
Variations due to the Virgo cluster and the Great Attractor are indicated. Notice 
that rates corresponding to about one event per year can be obtained only if the 
detector can probe distances at least of the order of 90-100 Mpc. 

The strength of a signal observed in a given detector is characterized by its signal- 
to-noise ratio S/N, which measures the signal amplitude in terms of the detector's 
noise characteristics. This quantity depends on the source power spectrum, the 
noise spectral density of the detector Sn{v) in Hz^^ and the adopted method of 
data analysis. For merging NS-NS pairs, the optimum S/N ratio is obtained by the 



D (Mpc) 



Fig. 3. Total expected NS-NS coalescence rate inside a volume of radius D as a function of D 



5, 2008 8:27 WSPC/INSTRUCTION FILE ijmp 



Expected coalescence rates of NS-NS binaries for laser beam interferometers 11 
SO called matched filter, namely, 



I ... (4) 



where | h{i') \^ is the sum of the square of the Fourier transform of both polarization 
components. According to Ref. 43, the equation above can be written as 

| = 8e(^)(^)^/6C(-™.x) (5) 

where D is the distance-luminosity to the source and the parameter tq is given by 
the relation 

ro = 9.25 X 10-^2 Mpc (6) 

with 

where Vq = 202.38 kHz. The term C(^'max) is defined by 

The spiral-in phase ends when the pair separation is such that tidal effects deform 
and may even disrupt the stars or when the last stable orbit is reached. In Eq. (8), 
i^max ~ 750 Hz corresponds to the maximum orbital frequency attained nearly the 
inner most circular orbit (see, for instance, Rcf. 44). The other parameter appearing 
in Eq. (5) is the "chirp" mass M = n^^^M^^^, with n and M being respectively 
the reduced and the total mass of the system. The angular function 6 depends on 
geometrical projection factors of the detector, F+ and Fx, for both polarization 
components (sec, for instance, Rcf. 45) and on the inclination between the orbital 
angular momentum of the binary and the line of sight, namely, 

=.4:[Fl{l + cos^if +4:F]iCos^i] (9) 

The probability distribution of Q was discussed in Ref. 46 and the authors have 
found an excellent approximation, which will be used here to estimate mean values 

P,e) = 5M!i,„<e<4 ,10) 

. /^From these equations, for a given S/N ratio and detector sensibility, the distance 
D probed by the instrument can be evaluated. For VIRGO and LIGO we have used 
the noise spectral density as given respectively by Refs. 47, 7 in a polynomial form, 
while for advanccd-LIGO we have adopted the expression given by Rcf. 48. Under 
these conditions, the parameter ro for each detector is 7.6 Mpc for VIRGO, 8.0 Mpc 
for LIGO and 120 Mpc for advanced-LIGO. Adopting S/N = 7.0, typical threshold 
for a false alarm rate of about one per year, the maximum probed distances are 13, 
14 and 207 Mpc for VIRGO, LIGO and advanced-LIGO respectively. Using data 
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from Fig. 2, the mean expected event rates are one each 125 yr for LIGO, one 
each 148 yr for VIRGO and 6 events per yr for advanced-LIGO. In a recent study, 
Refs. 49, 50 have considered the spectral noise density of VIRGO in the context of 
chirp signal detection. Four main spectral regions were identified, according to the 
dominant source of noise: v < 2 Hz (seismic noise); 2 Hz < v < 52 Hz (pendulum 
thermal noise); 52 Hz < i/ < 148 Hz (mirror thermal noise) and above 148 Hz (shot 
noise). Rather than directing efforts to the lowest frequencies in the attempt to 
observe a large number of cycles, these analyses indicate that the reduction in the 
mirror thermal noise band provides the highest gain in the S/N ratio. The direct 
dependence of the mirror thermal noise power spectrum on the temperature suggests 
the use of already existing cryogenic techniques. In parallel, also the pendulum 
thermal and shot noises should be reduced in the mirror thermal noise band. If 
the total noise in such band could be reduced by a factor of 10, a gain by a factor 
of 8 can be obtained, increasing the maximum probed distance up to 100 Mpc 
for S/N = 7 and allowing an expected detection rate of about 3 events every two 
years. A recently proposed sound configuration^^ for VIRGO would have a noise 
reduction through a large bandwidth and a similar expected detection rate. In this 
configuration the pendulum noise would reduced by a factor of 8, the mirror thermal 
noise by a factor of 7 and the shot noise by a factor of 4. These results suggest that 
only the planned sensibility of the second generation of interferometers will permit 
to probe the local universe deep enough in order that at least few events should 
be seen in a monitoring period of one or two years. An alternative strategy for 
the present generation of laser beam interferometers is the search of these events 
by using a network of detectors, since for a given false-alarm rate, the detection 
threshold is lowered as the number of detectors incrcascs^^^'''''. Data analysis of 
chirp signals detected by a network of detectors, in particular the one constituted 
by the two LIGO (Hanford and Livingston) and VIRGO, was recently considered by 
Ref. 54. They have estimated the increased detection sensitivity for this particular 
combination of detectors and for a false-alarm rate equal to one per year and a 
detection probability of 95%. Using the numbers of their Table HI, we derived a 
maximum probed distance of 22 Mpc and an expected detection rate of about one 
event every 26 yrs. 

4. Conclusions 

In spite of the numerous investigations performed by different groups in the past 
years on the coalescence rate of NS-NS binaries in the Galaxy, the uncertainties 
are still quite large, due to several badly known aspects concerning the evolution 
of massive binaries. In this work we have adopted a new approach to tackle this 
problem, calculating the coalescence rate from the star formation history of the 
galaxy and from the probability per unit of time for a NS-NS system to merger in a 
timescale r. In the case of the Galaxy, the star formation history was directly derived 
from observations, using the cromospheric activity index as an age estimator. The 
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data suggest that the star formation rate is intermittent, consistent with a picture 
in which the disk is gradually built by successive accretion episodes. Our numerical 
simulations on the evolution of massive binaries show that only 2.4% of the pairs 
remain bound after the second supernova explosion. Our simulations indicate also 
that the probability per unit of time for a NS-NS pair to coalesce in a timescale 
T is P{t) = B/t, in agreement with the early estimate by Ref. 20 and the recent 
simulations by Rcf. 36. Our approach indicate that the galactic coalescence rate is 
presently (1.7 ± 1.0) x 10^^ yr~^, a value which has remained more or less constant 
in the past 9 Gyr. Elliptical galaxies also contribute to the local mean coalescence 
rate. We have adopted as a typical E-galaxy a model from a grid able to reproduce 
the observed color-magnitude diagram of early-type galaxies in Virgo and Coma 
clusters. The star formation efficiency and the initial mass function, required to 
evaluate the fraction of formed massive stars, are fixed by fitting the photometric 
properties, whereas the fraction of formed massive binaries and the fraction of pairs 
which remain bound were taken to be equal to the galactic values. Under these 
conditions, the typical coalescence rate for an elliptical (or SO galaxy) of absolute 
magnitude Mb = -20.7 is about 8.6 x 10~^ yr~^. Notice that this value is about 
a factor of five higher than the galactic rate. The local coalescence rate, derived 
from an appropriate average weighted by the relative luminosities and abundances, 
is 3.4 X W~^yr~^. The expected coalescence within a volume of radius D was 
estimated by scaling the local mean rate with respect to the ratio between the 
total luminosity in the considered volume and the luminosity of the Galaxy. Galaxy 
counts from the LEDA data base were used to estimate the total luminosity. We 
have used the completeness estimates made by Ref. 38 and the Schechter luminosity 
function to correct for the absence of faint galaxies. The contribution of the Great 
Attractor was included by assuming that the total mass of this structure is about 
10^^ Mq and that the mean mass-to-luminosity ratio is 150 Mq/Lb,q- Our results 
indicate that a rate of about one event per year can only be obtained if detectors 
can probe distances at least of the order of 80 Mpc. This is not the case for the 
present generation of large laser beam interferometers, since from their planned 
sensibility, one event each 148 yr is expected for VIRGO whereas one event each 125 
yr is expected for LIGO. These rates arc modified when a network including three 
detectors (Hanford, Livingston and VIRGO) are considered, since in this case the 
expected rate corresponds to one event each 26 years. Recently proposed improved 
VIRGO configurations may significantly raise the expected detection rate up to 3 
events every two years. For advanced-LIGO, 6 events per year are expected to be 
seen. 
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